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HARMONIC  WAVE  SOLDTIONS  OF  IHE  NONLINEAR 
VORTICm  EQUATION  FOR  A ROTATING  VISCOUS  FLUID 

ABSTRACT 

With  the  assumption  of  non-divergent  horizontal  flow> 
harmonic  wave  solutions  of  the  nonlinear  hydrodynamic  equations 
for  a rotating,  viscous  fluid  are  obtained  both  for  the  plane 
and  for  the  sphere.  These  solutions  yield  waves  which  may  be 
damped,  amplified,  or  remain  unchanged,  depending  on  the  longi- 
tudinal and  latitudinal  extents  of  the  waves,  and  also  on  the 
vertical  profile  of  its  stream  function.  For  plane  flow,  the 
vertical  profile  of  the  velocity  components  is  shown  to  be  a 
sinusoidal  function,  or  a hyperbolic  sine  and  cosine  of  z, 
proportional  to  or  independent  of  z.  For  spherical  flow,  the 
vertical  profile  is  shown  to  be  proportional  to  a Bessel  func- 
tion of  order  n + A and  n,  a positive  integer,  is  the  order  of 
an  associated  Legendre  function.  Harmonic  waves  are  amplified, 
neutral  or  damped,  depending  respectively  on  whether  the  argu- 
ment of  the  Bessel  function  is  imaginary,  zero  or  real.  No 
stationary  wave  is  possible  in  spherical  viscous  flow,  although 
quasi-stationary  waves,  which  may  be  defined  as  waves  with  zero 
wave-velocity  but  changing  wave-amplitude,  may  occur  under  cer- 
tain circumstances. 


HARMONIC  WAVE  SOLUTIONS  OF  THE  NONLINEAR 
VORTICITY  EQUATION  FOR  A ROTATING  VISCOUS  FLUID 

1.  Introduction 

The  linearized  vorticity  equation  for  a rotating  invisid 
fluid  was  first  studied  by  Rossby  (1939),  who  obtained  a solu- 
tion for  the  motion  of  sinusoidal  waves  of  infinite  lateral 
extent  on  an  earth  whose  variation  of  Coriolis  paramter  with 
latitude  is  constant.  Later  Haurwitz  (1940)  extended  Rossby's 
method  and  obtained  solutions  for  waves  of  finite  lateral  ex- 
tent on  a horizontal  plane  (1940a)  and  on  a sphere  (1940b). 

Ertel  (1943)  treated  the  nonlinear  vorticity  equation 
for  a rotating  invisid  fluid  and  obtained  a solution  for  sta- 
tionary motion  of  harmonic  waves  of  finite  lateral  extent  on 
a sphere.  Craig  (1945),  Neamtan  (1946),  Rombakis  (1948), 
HBiland  (1951),  and  Long  (1952)  have  more  recently  studied 
the  nonlinear  vorticity  equation  in  great  detail. 

In  all  the  previous  treatments  the  fluid  has  been  consid- 
ered as  nonviscous,  and  motion  is  independent  of  height.  It 
is  the  objective  of  the  present  paper  to  study  the  fluid  motion 
including  the  effects  of  the  viscosity  and  the  height  depend- 
ency of  velocity  components.  On  account  of  frictional  effects, 
the  stream  functions  obtained  for  the  solutions  vary  with  time. 
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2.  Harmonic  Wave  Solutions  of  the  Nonlinear  Vorticity  Equation 
for  a (rotating)  Viscous  Fluid  on  an  Earth  whose  Variation  of 
Coriolis  Parameter  with  Latitude  is  Constant 

If  the  vertical  velocity  is  neglected  and  zero  horizontal 
divergence  is  assumed,  the  equations  of  atmospheric  motion  and 
equation  of  continuity  are 


diL 

at 


+ u 


8 u 
dx 


- (1) 


dx  ^ df  ' 


(3) 


where  x,  y,  z are  the  Cartesian  coordinates,  positive  toward 
the  east,  north  and  upward  respectively;  u,  v the  velocity 
components  in  the  x-  and  y-direction  respectively;  f = 2 A sin  ^ 
where  A is  the  angular  velocity  of  the  rotating  earth;  ^ is 
the  latitude;  v is  the  coefficient  of  kinematic  viscosity, 
p the  pressure  and  p the  density. 

Equation  (3)  shows  that  there  exists  a stream  function  >|i 
such  that 


di 


= ^ , 

8% 


(4) 


If  we  assume  further  that  fluid  is  autobarotropic,  the  vortlc- 


\ 


f 
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ity  equation  of  a rotating  viscous  fluid  is 


\at  a%  at  a^  ax  ^a7^  a^ 


where 


p._5S_cot*  , 


a being  the  radius  of  the  earth. 

Assume  wave  solutions  of  the  form 

C /co»  K(x-ct)co«Kijy\-»  (7a) 

♦ {x,y,i,t)=  -Uy  + ■<  A tin  (ay  + «)  + Bl  )>  • 

I \sin  K(x-v»-ct)  /J  (YIj) 

• g(z)e^'*  « 

(eotK(x-fRz->et)cotKt|y  (7c) 

cot  K(x-ct)  cot  K(i}y  + Rz)|\*(7d) 
tin  K(x-ijy«-Rz-ct)  ' 

•e-f.*, 


where  A,  B,  R,  < : arbitrary  constants 

(I,  R,  ' <^onstants  whose  values  are  to  be  determined 
c = wave  velocity,  positive  toward  the  east 
K = L = wave  length 

In  (7a)  (7c)  (7d)  : ’?  = k D = twice  the  lateral 

“ extent  of  the  waves 

In  (7b)  and  (7e)  : V = reciprocal  of  the  horizontal 

slope  of  the  trough  and  ridge 
lines 
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Stream  functions  (7a),  (7c),  and  (7d)  represent  waves  of  finite 
lateral  extent,  whereas  stream  functions  (7b)  and  (7e)  repre- 
sent waves  with  horizontally  tilted  troughs  and  ridges.  There 
is  vertical  tilt  of  troughs  and  ridges  in  the  waves  represented 
by  (7c),  (7d),  and  (7e)  but  no  vertical  tilt  in  the  waves  rep- 
resented by  (7a)  and  (7b). 

It  is  to  be  noticed  that  the  mean  zonal  speed  averaged 
over  a complete  wave  length, 

-L 

r udx  B U-«Acos  (ay<-«)  g(z)e~''^  , (8a,b) 

■'o 

dx  = U-«A coi(«y  + R*+f) , (8c,d,e) 

varies  not  only  with  latitude  but  also  with  height  and  time. 

Here  (8a),  (8b),  (8c),  (8d),  and  (8e)  give  the  mean  zonal  speed 
corresponding  to  the  waves  represented  by  (7a),  (7b),  (7c),  (7d), 
and  (7e),  respectively. 

Substituting  from  (7a),  (7b)  in  (5),  dividing  through 
with  g(z)e  and  arranging  terms,  we  have 


BK  I I^K^d  + ,*)  ( U - c)  + /ij  + Aa  - K^(l  + ,*)]  co*  (ay  + «) g (i) 


sinK(x-ct)cotKi}y 


co*K(x-ijy-c 


Acot(ay-^t) 


(9a) 

(9b) 


+ BK^I  + 1J*) 
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where  g"(z)  denotes  the  second  derivative  of  g(z).  Similarly, 
substitution  of  (7c),  (7d),  and  (7e)  in  (5)  gives 


1 1 + 1»*H  U - c)  + /s]  + A a [d*- 1 + 1J*)]  COI  («  y + R I + «)  •"  • 

.•lnK(x  + Ri-ct)  cot  Kijy. 

f*  n 

• f coiK(x-ct)  tin  K(ijy + Rz)  j + 

aH  f^-i»(d*+R*)  A tin  (ayt-Rz-fc) 

cotK(x-t}y  + Ri-ct) 

L>  J 

(9c) 

f*  *1 

.cot  K(x  + Rz-ct)cot  Kify  . 

+ BK*(I+/)  [e- vK^I+  1,%R^)J 

( cot  K(x-ct)cot  K(iyy+Rz)|  [ ■ 0 . 

(9d) 

\ i ^ 

tin  K(x-ijy  + Rz -ct) 

(9e) 

Thus  it  is  seen  that  (7a),  (7b),  (7c),  (7d),  and  (7e)  are  solu- 
tions of  (5)  providing  that 

. (10a,c,d) 

a » ■ 

From  (9a),  (9b),  (9c),  (9d),  and  (9e),  we  have 


2 2 

. .. 

u 2^ 

4»*'(L  +0^) 

(lla,c,d) 

.u--^L , 

4»  (i  + O 

(llb,e) 

(12a) 

(12b) 

(12c,d) 


V^[l*+d*(i  + r^] 

+ * (12e) 


and 

W = KR  (13) 

Equations  (11a,  c,  d)  give  wave  velocity  in  terms  of  U,  )8  , the 
wave  length  L,  and  the  lateral  extent  D/2,  whereas  (11b,  e)  give 
the  wave  velocity  in  terms  of  U,  )9  , L,  and  the  reciprocal  of 
the  horizontal  slope  of  trough  and  ridge  lines  of  the  waves. 
Since  , has  been  assumed  as  a constant,  g(z)  must  be  propor- 
tional to  the  hyperbolic  sine  and  hyperbolic  cosine  of  z,  a 
combination  of  sine  and  cosine  of  z,  proportional  to  z,  or  in- 
dependent of  z.  The  last  two  cases  contribute  no  frictional 
effect  and  are  of  little  interest.  We  shall  consider  the 
former  two  cases  which  will  be  shown  to  be  the  major  cause  of 
the  growth  of  waves  under  certain  circumstances.  Let 

g(z)  - G tinh  {x -I- E sinh  (z  , (14) 

where  G and  E are  arbitrary  constants.  It  is  obvious  g(z)  be- 
comes a circular  function  if  C is  imaginary. 

Stream  functions  (7a),  (7b),  (7c),  and  (7e)  may  respec- 


tively be  written 
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Uy  + 


r ■ /co»  K(x-ct)  coi  Kny\^ 

I A tin  (K^  1^4.  2 y-f«)-fB(  ||* 

'■  \$in  K(x -ijy -ct)  /•' 

/ 2 2 2 

•(GtinhCz  + EcoihCz)*'”'*^ 


(15a) 

(15b) 


f s - Uy  -(■ 


|a  tin(KA/i 


I + 1J*  y + KRz-t-«) 
cot  (x-f  Rz -ct)  cot  Ki;y 


■( 


(15c) 

\ 2 2 2 

+ B(  cot  K(x-ct)cot  K(ijy  + R»)jl  • (l  + ij+R  )t  ^ (15d) 
tin  K(x-ijy +Rz-ct)  (15e) 


The  first  term  of  the  right  side  of  equation  (15)  gives  a con- 
stant zonal  velocity,  the  second  term  a pure  zonal  current 
which  is  a function  of  latitude,  height  and  time,  the  last 
term  gives  both  zonal  and  meridional  velocity  which  is  a func- 
tion of  three  coordinates  and  time.  It  is  of  importance  to 
note  that  waves  will  always  be  damped  if  the  vertical  profile 
of  velocity  components  is  a sinusoidal  function  of  z,  propor- 
tional to  z,  or  independent  of  z.  If  the  vertical  profile  of 
the  horizontal  velocity  component  is  a combination  of  hyper- 
bolic sine  and  cosine  of  z,  waves  may  be  damped,  unchanged, 
or  amplified,  depending  respectively  on 

K“(l  + .*lj{*,  (16) 


\ 
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Examination  of  (16)  reveals  that  the  growth  of  waves  in  such 
a system  depends  on  the  balance  of  the  wave  length  and  the 
lateral  extent  (or  the  inclination  of  trough  lines)  of  the 
waves  with  the  vertical  profile  of  velocity  components.  The 
smaller  the  horizontal  extent  of  the  waves,  the  greater  the 
damping  effect  will  be.  The  frictional  effect  generated  by 
a sinusoidal  velocity  profile  always  serves  as  a damping  fac- 
tor, whereas  a hyperbolic  sine  and  cosine  profile  of  velocity 
serves  as  an  amplifying  factor.  This  can  be  justified  from 
the  following  energy  considerations.  Multiplying  (1)  (2)  by 
u and  V respectively,  then  adding  these  two  equations  and 
integrating  it  over  the  volume  of  the  system,  we  have 


Examining  the  above  equation,  one  finds  that  frictional  effect 
tends  to  decrease  the  total  kinetic  energy  of  the  wave  system 
if  the  velocity  components  of  the  system  are  sine  and  cosine 
functions  of  x,  y,  and  z;  whereas  it  tends  to  increase  the 
total  kinetic  energy  of  the  system  if  the  velocity  components 
are  hyperbolic  sine  and  cosine  functions  of  x,  y and  z.  ^ The 
^Friction  also  tends  to  increase  the  kinetic  energy  of  the 


\ 


velocity  components  of  the  solutions  with  which  we  are  deal- 
ing in  this  section  are  sine  and  cosine  functions  in  x and  y, 
but  may  be  independent  of  z,  proportional  to  z,  sine  and  co- 
sine functions  of  z,  or  hyperbolic  sine  and  cosine  functions 
of  z.  According  to  the  energy  considerations  it  is  obvious 
that  waves  will  always  be  damped  out  in  the  first  three  cases, 
since  the  total  kinetic  energy  of  the  system  will  decrease  due 
to  frictional  effect.  However,  in  the  last  case  the  change  of 
total  kinetic  energy  of  the  system  due  to  frictional  effect, 
and  therefore  the  growth  of  the  wave,  depends  on  the  balance 
of  the  kinetic  energy  being  fed  into  the  system  due  to  the 
vertical  hyperbolic  sine  and  cosine  profile  of  velocity  com- 
ponents, and  the  kinetic  energy  being  dissipated  due  to  the 
horizontal  sine  and  cosine  profile  of  the  velocity  components. 

It  has  been  previously  pointed  out  that  solutions  (15a, 
b,  c,  d,  e)  permit  not  only  vertical  but  also  meridional  vari- 
ations of  the  mean  zonal  current.  This  allows  us  to  approxi- 
mate a Jet  stream  by  a sinusoidal  vertical  profile,  which 
vanishes  at  the  earth's  surface  and  at  z = H,  say  the  height 
of  a homogeneous  atmosphere,  and  has  a maximum  speed  at 
z = H/2.  For  this  case  stream  functions  (15a)  and  (15b) 

system  if  u,  v are  proportional  to  x",  y",  z”  for  n5  2.  How- 
ever, this  case  is  not  of  interest  in  the  present  problem. 


become 
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♦ « |a  iln  (K-v/ l + ij*  » + «)  + B ( 

'•  '«tnK(x 


cot  K (x-et)eoi 
ijy-ct) 

A 

Z...  2. 


(18a) 

(18b) 


It  is  clear  that  in  such  a system  frictional  effects  tend  to 
damp  out  harmonic  waves.  Of  course  atmospheric  motion  is  com- 
plicated; wave  development  is  also  affected  by  pressure-density 
solenoids  and  vertical  motion.  Our  objective  here  is  to  bring 
out  the  effect  of  inertial  force  and  frictional  force  in  the 
atmospheric  motion. 

More  general  solutions  of  the  forms  (iSa,  b,  c,  d,  e)  may 
be  constructed  and  expressed  as  follows: 


( 


♦ ■-Uy+'(Ailn(dy + «|+ 


eo(  K|(x-ct)cot  Kjijjy 
tin  Kj(x-ij.y-ct) 


* (Gtinh  {z  -*-Ecosti 


Cz)e 


(19a) 

(19b) 


f » -Uy  +^|A.tin(tfy+K.Ri+«) 

i*r 

cot  Kj(x+Rz-ct)  cot 


-t-B, 


^cot  K|(x-ct)co»  Kj(ij|y  + Rz^J-( 
tin  Kj(x-i>|y+Rz  -et) 


2 2 2 
•»K,  (l+i},+R  )t 


(19c) 

(19d) 

(I9e) 


Equation  (20)  indicates  a connection  between  the  wave  length 
and  the  lateral  extent  in  (19a),  (19c)  and  (19d),  between  the 
wave  length  and  the  horizontal  inclination  of  trough-and-ridge- 
lines  in  (19b)  and  (l9e).  It  is  interesting  to  note  that  in 
viscous  plane  flow  stationary  waves  occur  if  both 

and  » 

£*= 

are  satisfied. 

Solutions  (19a),  (l9c)  and  (l9d)  contain  the  solutions 
of  Rossby  (1939),  Haurwitz  (1940a),  Craig  (1945),  and  Neantan 
(1946)  if  we  put  whereas  solutions  (19b)  and  (19e) 

contain  the  solutions  of  Machta  (1949)  and  Arakawa  (1953)  by 
putting  v3{sR«o. 

3.  Harmonic  Wave  Solution  of  the  Nonlinear  Vorticity  Equation 
for  a Rotating  Viscous  Fluid  on  a Sphere 

If  radial  velocity  is  neglected  and  zero  horizontal  di- 
vergence is  assumed,  the  two  fundamental  equations  of  spheri- 
cal motion  and  equation  of  continuity  are  respectively  reduced 


t 
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at  f as 


-p-  ton  9 


^ T ao  r as  T 2 


L r tin  9 r tin  6 o\  j 


|^(U,tin9)  + |^  = 0 , 


where 


yS  y-  ^ (,in  9-f^H-  y * r"fT  * 

7tin9  7tin9  « 5r 


Is  the  three  dimensional  spherical  Laplacian  operator.  d,X,r 
are  respectively  the  colatitude,  longitude,  and  radial  distance 
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from  the  center  of  the  earth;  0^,  are  respectively  the  ve- 
locity components  In  the  direction  of  the  gradient  of  6 and  X ; 
k and  H are  the  constants  of  universal  attraction  and  the  mass 
of  the  earth  respectively. 

Equation  of  continuity  (23)  shows  that  there  exists  a 
stream  function  ^ such  that 


u L_iJl 

• r tinfl 


Ux= 


r ae  * 


(25) 

(26) 


The  radial  component  of  relative  vorticity  may  be  written 


“ > It I.’ * -5^55 1?}  * • 

The  radial  component  of  absolute  vorticity  is 

q^=  2flcot#+ , 


where 


I 1 

tin^e  ai^j 


(27) 


(26) 

(29) 


is  the  spherical  surface-Laplacian  operator. 

If  we  further  assume  that  the  fluid  is  homogeneous^  the 
equations  of  motion  may  further  be  written  as  follows: 


\ 

\ 
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gu, 

at 


-Uj^(q+2  A cot  B) 


(30) 


where 


aui 

^^+U^(q+2A  cot0) 


(31) 


I a 

nine  5x 


[(u*+uJ)tA-M.XeViin'«jt»Hj, 


• '»  .Z.i.Zn  ax  ’ 


r Sin  9 r tin 


(32) 


r tin  fl  r tin  0 o'* 


(33) 


The  vorticity  equation  of  a rotating  viscous  fluid  on  a 
sphere  may  be  written 

2 

Ottf) 

(34) 


To  solve  (34)  we  assume  a wave  solution  of  the  form 
'Kfl.v.t)  * - wr^cot  0 ^l0M)f(r)  t'"**  , 


(35) 


where  w is  a constant,  Y^(0,X,t)  an  associated  Legendre  function 

of  order  n with  time  dependency,  s is  a parameter  which  will  be 

2 

shown  to  be  related  to  f(r).  The  term,  -ur  cos  6,  represents 


\ 
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the  part  of  the  stream  function  resulting  from  a solid  rotation 
of  angular  velocity  w relative  to  the  rotating  earth.  This 
term  is  independent  of  time  due  to  the  fact  that  solid  rota- 
tion is  not  affected  by  friction.  By  making  use  of  the  fact 

that  Y satisfies  the  equation 
n 

AY„+n(n+l)Y„»0  , (36) 

and  also  that: 

Aco<  Ss -2cotS  , (37) 

we  have  from  (35) : 

A^  + 20r*co«9  * 2tfl  + «)r%os#  — n(n-H)Y^f  (r)  • (38) 

After  substitution  of  (35)  and  (38)  into  (33),  the  vorticlty 
equation  reduces  to 


+ 


2(ft4-w)l 

.r,niff  ±U  -ilillll 

n(n+ljJ 

i5x 

f(r)Jj 

(39) 


It  is  evident  that  the  associated  Legendre  function  Yg(6,X,t) 
is  a solution  of  differential  equation  (39)  if 


and 


where 


Y„W,X,t)*  Y^W,X-ftnt)  , 


2(fl+i.>) 
n(n+l)  ' 


(40) 


(41) 

(42) 


I 
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Differential  equation  (41)  has  a solution 

iU)=rz  |(i*r)  , 

where  Z | is  a Bessel  function  of  order  n 1 • 

n+5  1 

Solution  (35)  is,  therefore 
♦ (fl,X,r,t)  = - itfr^cos  fl  + •^A|^P^(cos  fl)  s9)cos  m(X- 


•rZ  I 
n+-|f 

Where  Ajj  and  ^ (m  = 1,2,  ...  , n)  are  n + 1 amplitude  factors, 
o , 

(m  = 1,2,  ...  , n)  are  n phase  angles.  Whether  the  harmonic 
waves  of  this  system  are  amplified,  unchanged,  or  damped,  de- 


pends respectively  on 


S|0, 


10^) 


where  s is  associated  with  the  argument  of  Bessel  function 
4 X 

Z I (s*  r)  , For  large  values  of  r (in  the  atmosphere 

r = a + z and  s^o  is  of  the  order  of  10^)  '.r‘  Z i(s*r) 

’ n+Y 

behaves  like  a circular  function  or  an  exponential  function 
according  to  whether  s is  a positive  or  negative  number.  This 
is  similar  to  what  we  have  found  for  the  case  of  viscous  flow 


2 a'«6.37  X 10^  m,  s^~  ^-3  x lo^  . 

H 
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on  a plane.  It  is  also  interesting  to  note  that  s = 0 

(o+Z)^ 

if  f(r)  is  proportional  to  or  independent  of  r.  For  this  case 

flow  is  approximately  equivalent  to  non-viscous  flow. 

If  all  amplitude  factors  vanish  except  for  a certain 

m,  (45)  becomes 
2 

cot9 


f ) i ' 

+ < A„P„(co«fl)  + A^cotfl)  cot  m r Z i (•  r)i“  *'**, 

I*'*'  •'*'  nnj  (j  + i • 

for  m = 1,2)  ...  , n,  which  of  course  is  a solution  of  (34). 

It  is  seen  from  (46)  and  (42)  that  if  the  relative  angular 
velocity, «i , is  absent,  the  westward  wave  velocity  increases 
with  decreasing  wave  length  and  Increasing  lateral  extent  of 
the  waves.  The  corresponding  angular  velocity  of  the  current 
in  the  direction  of  the  gradient  of  8 and  X are,  respectively, 

dfl  , I S'* 

<•*  " d\ 


- m cos 


8^A'JJf^(eo$fl)sin  + 


Z 

ni- 


x' 

2 


(47) 


I 


X M 


dPn(cosfl) 

d(cos9) 


m dfflcosB) 
» d(cosB) 


cos  i+b!”) 
n n 


n-t- 


^ - 
X.U  r)« 


■vst 


(48) 


20  - 


Comparison  of  the  corresponding  terms  in  (46),  (47),  and  (48) 
shows  that  the  first  term  of  the  right  hand  side  of  (46)  is 
the  part  of  the  stream  function  resulting  from  a solid  rota- 
tion of  angular  velocity  w relative  to  the  rotating  earth. 

The  second  term  of  the  right  hand  side  of  (46)  gives  a purely 
zonal  current,  whereas  the  last  term  gives  both  zonal  and  me- 
ridional flows. 

The  part  of  the  stream  function  representing  the  last 
term  of  (46)  gives  2m  nodal  meridians,  or  m waves  along  each 
latitudinal  circle.  In  addition: 


P„(cos0)=sin  6 


dCcotS) 


This  is  a polynomial  in  cos  0 of  degree  (n  - ra)  which  vanishes 
on  n - m parallels  between  the  poles.  Therefore,  this  part  of 
the  stream  function  gives  cells  over  the  globe  which  is  total 
2m  (n  - m + 1)  in  number.  These  traveling  cells  possess  a 
period 


2» 


n(n  + l)  j 


m 


(49) 


It  is  interesting  to  note  that  invisid,  spherical-flow  waves 
become  stationary  if 


2ft 

(n-l)(n+2)  • 


W - - 


(50) 
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It  is  seen  from  (46)  that  no  stationary  wave  is  possible  in 

viscous  flow  over  a rotating  sphere,  since  this  condition  would 

also  require  s = 0.  But  Z | (0)=  0 for  n -1/2.  However, 

n+-^ 

quasi-stationary  waves,  which  may  be  defined  as  waves  with  zero 
wave-velocity  but  changing  wave-amplitude,  may  occur  when  con- 
dition (50)  is  satisfied.  In  order  to  illustrate  the  distri- 
bution of  harmonic  waves  over  a spherical  surface,  the  stream 
function  of  an  example  for  quasi-stationary  case  with  m = 4 
and  n = 5 has  been  sketched  as  shown  in  fig.  1. 

4.  Conclusions 

It  has  been  shown  that  solutions  of  the  nonlinear  vorti- 
city  equation  for  a rotating  viscous  fluid  yield  waves  which 
vary  with  both  height  and  time.  Whether  these  waves  are  damped, 
neutral,  or  amplified  depends  on  the  longitudinal  and  latitu- 
dinal extents  of  the  waves,  and  also  on  the  vertical  profile 
of  their  stream  function  (therefore  their  velocity  components). 
This  is  evidently  due  to  the  effect  of  Navier-Stokes  friction. 

Plane-flow  waves,  for  which  velocity  components  are  pro- 
portional to  or  independent  of  height  or  have  a sinusoidal 
vertical  profile,  appear  always  to  be  damped;  on  the  other 
hand,  waves,  for  which  velocity  components  vary  with  the  hyper- 
bolic sine  and  cosine  of  z,  may  be  amplified,  neutral,  or 


i 
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2 2 < 2 

damped,  depending  respectively  on  K (l  + ij  )?C 
locity  of  propagation  of  these  waves  is 


The  ve- 


c = u--^r:  , 

Kll+i;  ) 

which  is  independent  of  height.  However,  it  is  of  importance 
to  note  that  here  U is  not  generally  equal  to  the  mean  zonal 
velocity. 

For  spherical  flow,  waves  are  amplified,  neutral  or  damp- 
ed, depending  respectively  on  s = 0,  where  s is  associated 
with  the  vertical  profile  of  the  stream  function,  which  was 

i -k 

found  to  be  proportional  to  r^Z  i (s  r)  . For  large  values 

i i i 

of  s r,  r Z I (s  r)  behaves  like  a sinusoidal  or  exponential 
function  of  s^r  , depending  respectively  on  whether  s is  a 
positive  or  negative  quantity.  The  period  of  these  spherical 


waves  was  found  to  be 


r = — 

[ n(n  + l) 


It  was  also  found  that  no  stationary  wave  is  possible  in 
spherical  viscous  flow,  although  quasi-stationary  waves,  which 
may  be  defined  as  waves  with  zero  wave-velocity  but  changing 
wave-amplitude,  may  occur  when 


(n-l)(n+2)  * 


f 


"ar 


^ \ 
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Harmonic  waves  are  generally  damped  in  a laminar  viscous 
flow  between  two  rotating  plates  or  spherical  shells,  since 
waves  must  vanish  on  the  solid  boundaries.  This  permits  only 
a sinusoidal  vertical  profile  of  the  stream  function.  It  is 
seen  from  the  properties  of  Bessel  function  that  for  a certain 
integer  n the  greater  the  value  of  $,  the  smaller  the  distance 
between  neighboring  zeros  of  Z , (s^r)  . This  means  that 

the  smaller  the  distance  between  two  neighboring  zeros  of  the 
vertical  profile  of  the  stream  function,  the  greater  is  the 
damping  factor. 
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